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Abstract 

Order statistics find applications in various areas of communications and signal processing. In this 
paper, we introduce an unified analytical framework to determine the joint statistics of partial sums of 
ordered random variables (RVs). With the proposed approach, we can systematically derive the joint 
statistics of any partial sums of ordered statistics, in terms of the moment generating function (MGF) 
and the probability density function (PDF). Our MGF-based approach applies not only when all the K 
ordered RVs are involved but also when only the K s (K s < K) best RVs are considered. In addition, 
we present the closed-form expressions for the exponential RV special case. These results apply to the 
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performance analysis of various wireless communication systems over fading channels. 

Index Terms 

Joint PDF, Moment generating function (MGF), Order statistics, Probability density function (PDF), 
Rayleigh fading. 

I. Introduction 

The subject of order statistics deals with the properties and distributions of the ordered 
random variables (RVs) and their functions. It has found applications in many areas of statistical 
theory and practice [1], with examples including life-testing, quality control, signal and image 
processing [2], [3]. Recently, order statistics makes a growing number of appearance in the 
analysis and design of wireless communication systems (see for example [4]-[14]). For example, 
diversity techniques have been used for over the past fifty years to mitigate the effects of fading on 
wireless communication systems. These techniques improve the performance of wireless systems 
over fading channels by generating and combining multiple replicas of the same information 
bearing signal at the receiver. The analysis of low-complexity selection combining schemes, 
which select the best replica, requires some basic results of order statistics, i.e. the distribution 
functions of the largest one among several random variables. 

More recently, the design and analysis of adaptive diversity combining techniques and mul- 
tiuser scheduling strategies call for some further results on order statistics [5], [6]. In particular, 
the joint statistics of partial sums of ordered RVs are often necessary for the accurate character- 
ization of system performance [7], [14]. The major difficulty in obtaining the statistics of partial 
sums of ordered RVs resides in the fact that even if the original unordered RVs are independently 
distributed, their ordered versions are necessarily dependent due to the inequality relations 
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among them. Recently, the co-author has applied a successive conditioning approach to convert 
dependent ordered random variables to independent unordered ones [5], [6]. That approach, 
however, requires some case-specific manipulations, which may not always be generalizable. 

In this paper, we present an unified analytical framework to determine the joint statistics of 
partial sums of ordered RVs using a moment generating functions (MGF) based approach. More 
specifically, we extend the result in [15]— [17], which only derive the joint MGF of the selected 
individual order statistics and the sum of the remaining ones, and systematically solve for the 
joint statistics of arbitrary partial sums of ordered RVs. The main advantage of the proposed 
MGF-based unified framework is that it applies not only to the cases when all the K ordered 
RVs are considered but also to those cases when only the K s (K s < K) best RVs are involved. 
After considering several illustrative examples, we focus on the exponential RV special case and 
derive the closed-form expressions of the joint statistics. These statistical results can apply to 
the performance analysis of various wireless communication systems over generalized fading 
channels. 

The remainder of this paper is organized as follows. In section II, we summarize the main idea 
behind the proposed unified analytical framework, including the general idea and some special 
considerations. We then introduce some common functions and useful relations in section III, 
which will help make the results in later sections more compact. In section IV and V, we 
present some selected examples on the derivation of joint PDF based on our proposed approach. 
Following this, we show in section VI some closed form expressions for the selected examples 
presented in previous sections under i.i.d. Rayleigh fading conditions. Finally, we discuss some 
useful applications of these results in section VII. 
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II. The Main Idea 

Let oo > 7i:A- > 7 2: ftr > 73 : x • ■ ■ > ^k-.k > be the order statistics obtained by arranging 
K nonnegative i.i.d. RVs, {ji} i=1 , in decreasing order of magnitude. The objective is to derive 
the joint PDF of their partial sums involving either all K or the first K S (K S < K) ordered RVs, 

K m K s 

e.g. the joint PDF of j m:K and J2 7kk or the joint PDF of £ j n:K and £ In-.K- 

n=l n=l n=m+l 

A. General steps 

The proposed analytical framework adopts a general two-step approach: i) obtain the analytical 
expressions of the joint MGF of partial sums (not necessarily the partial sums of interest as will be 
seen later), ii) apply inverse Laplace transform to derive the joint PDF of partial sums (additional 
integration may be required to obtain the desired joint PDF). To facilitate the inverse Laplace 
transform calculation, the joint MGF from step i) should be made as compact as possible. An 
observation made in [15]— [17] involving the interchange of multiple integrals of ordered RVs 
becomes useful in the following analysis. Suppose for example that we need to evaluate a multiple 
integral over the range 7 a > 7i > 72 > 73 > 74 > 7&- More specifically, let 

rla rfi rj2 rls 

I= d 7i / d J2 / d-f 3 / dj A p (71, 72, 73, 74)- (1) 

J-f b J 7b J 7b Jib 

It can be shown by interchanging the order of integration, while ensuring each pair of limits is 
chosen to be as tight as possible, the multiple integral in (OQ) can be rewritten into the following 
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equivalent representations, 

rla rla rla rla 

1=1 dj4 dj 3 / d-f 2 I dji p(7i>72,73,74) 

J lb J 14 J 13 J 11 

rla rli rl3 rla 

= dj 2 dj 3 / d-f 4 / d7i p (71,72, 73) 74) 

Jib Jib Jib J 11 

rla rla rl3 rli 

= dj 3 dji / dj 4 / dj 2 P (71,72,73, 74) 

Jib J 13 Jib Jl3 

rla rli rli rli 

= / CZ71 / ^74 / dj 3 / c?72 p (71, 72, 73, 74) 

7b *' 7b J 14 J 13 

rla rla rli rli 

= dj 4 dn/x j dj 3 / dj 2 P (71, 72, 73, 74)- (2) 
J~f b J74 J 74 J73 

The general rule is that the integration limits should be selected as tight as possible using the 

remaining variables. For example, in the first equation of ©, the variables are integrated in the 

order of 71, 72, 73, and 74. Based on the given inequality condition 7„ > 71 > 72 > 73 > 74 > lb, 

the integration limit of 71 should be from 72 to 7 a because 72 is the tightest among the remaining 

RVs. Similarly, the integration limit 73 is from 74 to 7 a because 71 and 72 were already integrated 

out. 

After obtaining the joint MGF in a compact form, we can derive joint PDF of selected 
partial sum through inverse Laplace transform. For most cases of our interest, the joint MGF 
involves basic functions, for which the inverse Laplace transform can be calculated analytically. 
In the worst case, we may rely on the Bromwich contour integral. In most of the case, the 
result involves a single one-dimensional contour integration, which can be easily and accurately 
evaluated numerically with the help of integral tables [18], [19] or using standard mathematical 
packages such as Mathematica and Matlab. 

B. Special cases 

The general steps can be directly applied when all K ordered RVs are considered and the RVs 
in the partial sums are continuous. When these conditions do not hold, we need to apply some 
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extra steps in the analysis in order to obtain a valid joint MGF. Specifically, when only the best 
K s {K s < K) ordered RVs are involved in the partial sums, we should consider the J\ s th order 
statistics ^k s -.k separately. Without this separation, we cannot find the valid integration limit in 
calculating the joint MGF. As an illustration, the example in Fig. \T\ considers 3-dimensional joint 
PDF of the partial sums of three group of RVs which are {'fx-.K, Ii-.k, Iz-.k}, {ia-.k, 75:K, Iq-.k}, 
and {77:x, Js-.k}, with K = 10 and K s = 8. Following the proposed approach, we will derive 
the 4-dimensional joint MGF in step i) by considering 7 8:i< - separately, i.e. the joint MGF of 
the following four groups of RVs {^i :K , 12-.K, Is-.k}, {1a-.k, Js-.k, 1g-.k}, {1i-.k}, and {^s-.k}- 
After obtaining the corresponding 4-dimensional joint MGF, we need to perform another finite 
integration to obtain the desired 3-dimensional joint PDF. 

When the RVs involved in one partial sum is not continuous, i.e. separated by the other RVs, 
we need to divide these RVs into smaller sums. The example in Fig. [2] illustrate this process. If 
we consider 3-dimensional joint PDF of {^i :K , 12-.K, Is-.k, 1&-.k}, {iz-.k, Ia-.k}, and {'Jt-k^s-.k}, 
then there are three partial sum of RVs {j hK , j 2 . K , j 5 . K , 7 6:X }, {>y 3:K , j A:K }, and {-Jj-.k^Is-.k}- 
Note that the first group is split by the second group of RVs (as such discontinuous). More 
specifically, the second group, {73:^,74:^}, split the original group, {71^, ^-.k, Iz-.k, Ig-.k}, 
into two groups as {'Ji-.k , J2-.k}, and {75^, 'Jg-.k}- Therefore, we also consider this split group 
{7i:X,72:^} and {75^, 76 : x} as two smaller groups. As a result, we will derive 5-dimensional 
joint MGF of {-y 1:K , j 2 -.k}, {l3-.K, Ia-.k}, {Js-.k, Jg-.k}, {j7-.i<}, and {7 8: a'} in step i). Similarly 
to the first example, after the joint PDF of the new substituted partial sums are derived with 
inverse Laplace transform in step ii), we can transform it to a lower dimensional desired joint 
PDF with finite integration. 

The proposed approach is summarized in the flowchart given in Fig. [3j where we consider 
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three different case separately. In the following sections, we present several examples to illustrate 
the proposed analytical framework. Our focus is on how to obtain a compact expression of the 
joint MGFs, which can be greatly simplified with the application of the following function and 
relations. 

III. Common Functions and Useful Relations 

In this section, we introduce some common functions and their properties. These results will 
be used to simplify the derivation of joint MGFs in later sections. 

A. Common Functions 

i) A mixture of a CDF and an MGF c (7, A): 



where p (x) denotes the PDF of the RV of interest. Note that c (7, 0) = c (7) is the CDF 
and c(oo, A) leads to the MGF. Here, the variable 7 is real, while A can be complex, 
ii) A mixture of an exceedance distribution function (EDF) and an MGF, e (7, A): 




(3) 




7 



(4) 



Note that e (7, 0) = e (7) is the EDF while e (0, A) gives the MGF. 



iii) An Interval MGF /i {^ a , j b , A): 




(5) 



Note that fj, (0, 00, A) gives the MGF. 
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Note that the functions defined in ©, © and © are related as follows 

c( 7 ,A) = e(0,A)-e( 7 ,A) (6) 

= c (oo, A) — e (7, A) 
e( 7l A) = c(oo,A)-c(7,A) (7) 

= e (0, A) - c (7, A) 
jU(7 ,76,A) = c(76,A) -c(7a,A) (8) 

= e(7 Q ,A) -e(7b,A). 

B. Simplifying Relationship 

i) Integral I m : 

Based on the derivation given Appendix UH the integral I m defined as: 

Jo Jo 

X/ d'y m+ 2:KP(l,n+2:K)exp(X'f m+ 2:K) ■ ■ ■ d^ K :KP {l K-.k) {^-1 K:k) , (9) 

Jo JO 

can be expressed in terms of the function c (7, A) as 

/m = ,K \^ ^-l-.K, A)] (A '- l+1) . (10) 

(a — m + 1)\ 

ii) Integral I' m : 

Following the similar derivation as given in Appendix HH the integral I' m defined as 

pOO nOO 

I'm = / d-f m:K p (-f m:K ) CXp (X-f m:K ) dj m - 1:K p(j m -l:K)cxp{\jm-l:K) 

Jfm-\-l:K J*1m:K 

j* OO />oo 

X/ dj m -2:KPhm~2:K)cX-p(Xj m ^2:K) ■ ■ ■ / d-f 1:K p (j 1:K ) CXp ( A7l : if ) , (11) 

can be expressed in terms of the function e (7, A) as 



I'm = —j [e (7m+l:if, A)]™ . (12) 
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iii) Integral 7" 6 : 

Based on the derivation given in Appendix HTJl the integral I" b defined as 

I'a,b= \ d~f b -l:K p{lb-l:K)cKv(\~1 b -l; K ) I d^ b -2:KP{lb-2:K)cxp{\^b-2:K) 
Jlb-.K Jlb-l-.K 

fla-.K rla-.K 



X/ dj b - 3 ; K p (j b - 3 .. K ) CXp (\j b ^ 3:K ) ■ ■ ■ / dj a +l:KPha+l:K)cxp(Xj a+1 .. K ), (13) 

•>Jb-2:K J-1a + 2:K 

can be expressed in terms of the function /j (•, •) as 

1 

(b-a-l)l 



I 'U = 7T 1 -^T7lv(jb:K,Ja:KA)f ^ for 6 > a. (14) 



IV. Sample Cases when All K Ordered RVs are Considered 

Assume the original RVs {7^} are i.i.d. with a common arbitrary PDF p (7), the A'-dimensional 
joint PDF of {ji-.K}f = i is simply given by [1] 

K 

P(jl-.K,J2;K, ■ ' " ,1K-.k) = F ■ JJp (ji:K) for Jl:K > 72-.K > l3:K ■ > JK-.K, (15) 

i=l 

where F = K\. 

K 

Theorem 4.1: (PDF of ^ n , K among K ordered RVs) 

71=1 

K 

Let Zi = ln;K for convenience. We can derive the PDF of Z = [Zi] as 

71=1 

Pz ( Zl ) = C^{[c( o,-S 1 )] K }, (16) 

where C$1 {■} denotes the inverse Laplace transform with respect to Si. 
Proof: The MGF of Z is given by the expectation 

OO "11: K 

MGFz (Ai) = E{cxp(AiZi)} = F -JdjuKP (7i:/c) ex P (^i7v.k) ^Jdry^KP (12-.k) exp (Ai7 2: jr ) 

o o 

JK-1.K 

x-- - xJd>y K:K p('y K:K )exp(\i"/ K .. K ), (17) 
b 

where E{-} denotes the expectation operator. By applying (flOl) . we can obtain the MGF of 

Z l = Em=l lrn:K aS 

MGFz (Ai) - [cCoo.Ai)]*. (18) 
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Therefore, we can derive the PDF of Z\ = J2m=i Im-.K by applying the inverse Laplace transform 
as 



Pz { Zl ) = c£ {MGF Z (Si)} 



C s l i {[c{^-S 1 )] K ). (19) 



K 

Theorem 4.2: (Joint PDF of j m -.K and Iu-.k) 

n=l 

K 

Let Zi = 'jm-.K and Z 2 = Yl In-.K f° r convenience. We can obtain the 2-dimensional joint PDF 

n=l 

of Z = [Z U Z 2 ] as 

Pz(z\,z 2 ) =p k (zi,z 2 ) 

7m:K) 2~2 1n:K 
n=l 



£^yp(z{)£ s * |[c(zi,-S , 2 )] (A ' 15 1 for m = 1, zi > ^jz 2 

w -^— W p{ Zl )Cs 1 2 {[c(z x -S 2 )f K -^ [6(^,-52)]^} form > 2. 



(20) 



Proof: The second order MGF of Z = [Z±, Z 2 ] is given by the expectation 

MGF z (\ 1 ,\ 2 )=E{exp(\ 1 Z 1 +\ 2 Z 2 )} 
= F- 



j djuKP (Ii-.k) exp (A 2 7l:Jf ) jdj 2:K P (j 2: k) exp {X 2 J 2:K ) • • 7^ rf7m-l:KP (7m-l:Jf) CX P (A27 


7m-l:if 

X / d~1 m :KP (im-.K) exp (A 1 7 m: jf) 



m-l:K) 





7m:Jf 7K-1:K 

^m+iiA-p (7m+i:if ) exp (A 2 7 m +i : jir) • • ■ / d-f K:K p (j K . K ) exp (A 2 7,fi::_ftT). (21) 



We show in Appendix [IV] that by applying (flOl) . © and (fl"2l) . we can obtain the second order 

K 

MGF of Z x = Im-.K and Z 2 = 1™k as 

n=l 



MGF Z (X 1 ,X 2 ) 
F 



(K — m)\ (m — 1)! 



d-jm-.KP (im-.x) exp (Ai7„ 1: k) [c (7 m :if, X 2 )}^ K m ' [e (jm-.K t \ 2 )} {m ~ 1) . (22) 
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With the MGF expression given in (1221) at hand, we are now in the position to derive the 

K 

2-dimensional joint PDF of Z\ = j m: x and Z 2 = In-.K- Letting Ai = —S\ and A 2 = — S 2 , 

n=l 

we can derive the 2-dimensional joint PDF by applying the inverse Laplace transform as 



Pz ( Zl ,z 2 ) = Csl A {MGF Z {-Si, -S 2 )} 
K\ 



(K — m)\ (m - 1)! 



V {lm:K) C> S \ {exp (-5l7m:if )} 



XC 



S2 X {[c(l m :K,~S 2 )] {K - m) [e( 7m:K ,-S 2 )}^} 



Based on the inverse Laplace transform properties in Appendix HI 



(23) 



{exp {-Sij m:K )} = 6{zi- jm-.x) 



(24) 



Therefore, substituting (|24|) in (1231) we can obtain the desired 2-dimensional joint PDF of Z\ 



K 

j m:K and Z 2 = ln:K- ■ 

n=l 

m K 

Theorem 4.3: (Joint PDF of In-.K and Yl In-.K) 

n=l n=m+l 
rn K 

Let Z\ = Yl In-.K and Z 2 = In-.K for convenience, then we can derive the 2-dimensional 

n=l n=m+l 

joint PDF of Z = [Z U Z 2 ] as 



PZ (Zl, Z 2 ) = P m K (Z1,Z 2 ) 

n = l n = m+l 

= Cg^ Sa {MGF Z (—Si, —S 2 )} 

oo 

K\ 



{K-m)\(m-l)\ 



K 



l(jrn:K)£ S i { CX P (~SlJm:I<) [e {lm:K , "Si )] ( '" ^ } 



xCs 2 1 {[c( lm ;K,S 2 )}^- m) } 



for z\ > 



K — m 



z 2 . 



Proof: The second order MGF of Z = [Z±, Z 2 ] is given by the expectation 



(25) 



MGF Z (Xi,X 2 ) = E {exp (A t Zi + X 2 Z 2 )} 

= Fj drfuKP (iuk) exp (\ijuk) • • • J dj m:K p (7m:*:) exp (Xij m:K ) 



7m:K 1K-1-.K 

X /" d7m+l:J<rP (7m+l:if ) exp (A 2 7m+l : Rr) • • • J djK:KP(lK:K)exp(X 2 jK:K)- 




(26) 



(27) 
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We show in Appendix IVl that by applying (flOl) and © and then (fl"2~l) . we can obtain the second 
order MGF of Z as 

MGF Z (Ai,A 2 ) 

oo 

Ai)]^- 1 ). (28) 



(K - m)\ (m - 1)! 



Again, letting Ai = — S\ and A2 = —52, we can obtain the desired 2-dimensional joint PDF 

m K 

of Zi = ^ j n:K and Z 2 = Yl 7n-.K by applying the inverse Laplace transform. ■ 

n=l n=m+l 

V. Sample cases when only K s ordered RVs are considered 

Let us now consider the cases where only the best K s {< K) ordered RVs are involved. 
Assuming the original {7^} are i.i.d. RVs with a common arbitrary PDF p (7) and CDF P(j), 
the K s -dimensional joint PDF of {juk}^ is simply given by [1] 

P{jl-.K,J2:K, - ■ ■ ,JK s :K) = F - Y[p(-fi:K)[P {■yK s :K)} K ~ Ka , (29) 

i=l 

where F = K.\{*) = 

K s 

Theorem 5.1: (PDF of £ In-.K, K s > 2) 

71=1 

K s 

Let Z' = ^2 In-.K for convenience, then we can derive the PDF of Z' as 



71=1 



Pz> (x) = Pj2^i 1 7„. K ( x ) = J Pz{x- z 2 ,z 2 ) dz 2 for K s > 2, (30) 

where 

Pz *a) = {Ks F _ p (z 2 ) [c (z 2 )] ( " s) C~ s l {[e (z 2 , -S 1 )] {K ^} . (31) 

K„-l 

Proof: We only need to consider ^k s -.k separately in this case. Let Z 1 = ^ n . K and 

71=1 

Z2 = 1k 3 :K- The target second order MGF of Z = [Zx, Z 2 ] is given by the expectation in 

MGF z (Xi.X 2 ) = E{exp(AiZ! + \ 2 z 2 )} 

= F I d'yi-.xP (ii-.k) exp (Ai/yi^) • • • / dj Ks -i-.KP (jk s -i-.k) cxp (Xi^k.-i-.k) 

Jo Jo 

K S —1:K 



X 



d'fK s :KPhK s :K)exp(X 2 '-f Ks:K ) [c("ik s :k)] K s ■ (32) 
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By simply applying © and then (fl"2l) to (|32|) . we can obtain the second order MGF result as 



F 



MGF Z (Ai,A 2 ) = (v - I djK s :KP(7K s :K)exp(\ 2 jK s :K)[c{-fK s :K)} K K " [e{.lK 3 :K,M)] Ks 1 • (33) 

o 



Again, letting Ai = —S\ and A 2 = — S 2 , we can obtain the 2-dimensional joint PDF of Z\ = 

K s -1 

Jn-.K and Z 2 = Jk s :K by applying the inverse Laplace transform. Finally, noting that Z' = 

71=1 

Zi + Z 2 , we can obtain the target PDF of Z' with the following finite integration 

Pz'(x) = pz (x - z 2l z 2 ) dz 2 . (34) 



i) 



Theorem 5.2: (Joint PDF of •y m :K and ^u-.k) 

n=l 

Let X = 'jm-.K and Y = Yl ln-.K, then the joint PDF of Z = [X, Y] can be obtained as 

71=1 



Pz(x,y)=p k s (x,y) 

1m:K, J2 ln:K 
71=1 



S(K a s-2\ X p Ks-i (x, zz, y - z 2 ) dz 2 , m = 1, 



n = 2 



=< 



Jo I(m-i)x m)Z4 P" 1 - 1 K ° - 1 (z 1 ,x,y-z 1 -z 4 , za) dz x dz^ 1 < m < K s - 1, 

£ 7n:K,lm:K, ln:K ,7 K S :K 

n — 1 n—m + 1 (35) 

/r?f _2) as P*-- a (zx,x,y- zi)dz lt m = K a -l, 

Z) 7»:K:7jf s -l:Jf,7iST s :K 

71=1 

p jc,-i m = K s , 

lKs:K, 7m:X 
n=l 
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or equivalently 



P k s (x,y) 

lm:K, ln:K 

n=l 
n=im 



r{K s -2)x , \ , 

.\(k s -2\ P Ks-i {x,z 2 ,y- z 2 )dz 2 , 



n = 2 



Jo I(K s S -m-l)z 4 K 4^ 1 {y-Z 3 -Z4, X, Z 3 , Z 4 ) dz 3 dz 4 , 1 < 171 < K s - 1, 



5Z 7»:Jr>7m:K, E ln:K,lK s :K 
n—l n — m+1 



(36) 



JoPks-2 (y - z 3 ,x,z 3 )dz 3 , m = K s -l, 

12 ImKflKs-l-.K-TfKsiK 



71=1 



p m = K s . 

lK a :K; H 7n:if 
>■ n=l 

Proof: To derive the joint PDF of 7 m: /< and 7„ : /<, we need to consider four cases i) 

n=l 

m — 1, ii) 1 < m < A' s — 1, iii) m = K s — 1 and iv) m = K s separately based on our unified 
framework. While for case iv), we can start with the second order MGF of 7 m: x and ^ n ., K 

n=l 

directly, we should consider substituted groups instead of the original groups for cases i), ii), 
and iii). More specifically, for cases i) and iii), we need to consider ^k s -.k separately as shown 



in Fig. |4(a)| and |4(c)| whereas, for case ii), as one of original groups is split by r y m: K, we should 
consider substituted groups for the split group instead of original groups as shown in Fig. |4(b)| 
As a result, we will start by obtaining a four order MGF for case ii) and a three order MGF for 
case i) and case iii). In all these cases, the higher dimensional joint PDF can then be used to 
find the desired 2-dimensional joint PDF of interest by transformation. 

Applying the results in ©, (flOl ), (fl"2l) and (fT4"l) . we derive in Appendix [VTI the following joint 
MGF for different cases 
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a. m — 1 

Ks-l 

Let Z x = 7i :A -, Z 2 = In-.K, and Z 3 = j Ks:K , then 

n=2 

MGF^(Ai,A2,A 3 ) 



-F / d~i Ka .. K p(~iK 3 .. K )eyLv(\^ Ka .. K ) [c{^ Ks :k)] {K Ks) 



x J d'yuKP (ji-.k) exp (Ai7i : g) ^ ^ ^ [yu ^jr,:*-, 71:1c, A 2 )] ( ^ a 2) . 
b. 1 < m < K s — 1 

m-l K B -\ 

Let Zi = XI 7«:/<' ^2 = Im-.K, ^3 = X 7n:^, and Z 4 = 7^ s:A , then 

n=l n=m+l 



MGF Z (Ai, A2, A3, A 4 ) 
F 

~ (K s -m - l)!(m- 1)! 



d^K s :KP {ik s -.k) exp (A47if s :x) [c (7nr,:ir)] 



(K-K 3 



x / rf7m:xp(7m:K)exp(A 2 7 m: K) [e(7 m:K ,Ai)] (m :) [M(7K s :/c,7m:K, A 3 )] (Ks m J) . 



c. m = K s — 1 

K a -2 

Let Zx= J2 -y n:K , Z 2 = Jk s -1:K and Z 3 = 7^^, then 

n=l 

MGFz (Ai, A2, A3) 

OO 

= F / rf7K- s :KP(7K s :K) G Xp(A 3 7K s :K) [c(7K s: K)] (i:f ~ i:fs) 



x y d7A- 3 -i:A-p(7A: s -i:K-)cxp(A27K- 3 -i:A-) jj£ 1 -2)\ ^ e ^ K ^- 1:K ' Xl ^ Ks 2) ■ 

d. m = K s 

K s -1 

Let Zi = 7^ s :/< and Z 2 = |] 7^, then 

n=l 

00 

MGFz (Ai, A2) = F J d-yK s :KP{lK s :K)exp(^i7K s :K)[c(iK 3 :K)] (K ^ Ks) 





X 



(A s - 1)! 



Starting from the MGF expressions given above, we apply inverse Laplace transforms in 
pendix [Vj in order to derive the following joint PDFs 
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a. m — 1 

Pz (zi, z 2 ,z 3 ) = p Ks -i (z 1 ,z 2 ,z 3 ) 

n = 2 

for z 3 < zi, (if, - 2) z 3 < z 2 < (K s -2)zx. (41) 



where U (•) is the unit step function, 
b. 1 < m < K s — 1 



PZ {Zl, Z2, Z3, Z&) = Pm-i K,-i (z 1 ,Z 2 ,Z 3 ,Z 4 ) 

22 "fn-.Kifm.K, Yl ln:K,lK s :K 

n.=l n=m+l 

P(Z 2 )P(^4) [c(z 4 )] (X ~ Ks) t/(z 2 - Z 4 ) 



(isT s -m - l)!(m- 1)! 
x£s 1 1 {[e(^,-5i)] (w - 1) }^, 1 {[^(^^-S3)] (Jf '- ro - 1) }, 
for Z4 < Z2j z\ > (m — 1)^2 and (if s — m — 1) 2:4 < Z3 < (K s — m — 1) z 2 - (42) 

C. m = i^. s — 1 

(^1,^2,^3) =Pk s -2 (zi,z 2 ,z 3 ) 

Yl TmK'jYK's — 1:Kj7K s :K 



F 



(z 2 )p(z 3 ) [c{z 3 )] [K - K ^ U (z 2 - z 3 ) Csl {[e (z 2 , -S^*' -35 } , 



(K,-2)! J 

for z 3 < z 2 , 21 > (A" s - 2) z 2 . (43) 



d. m = 



Pz{zi,z 2 )=p k.-i (zi,z 2 ) 

1K S .K, Y ln:K 



F 



- 1)!' 

for 2 2 > (44) 
Finally, the joint PDF of j m .K an d Yl In-.K can be obtained. ■ 

n=l 

Note that (1351) and (l36l) involve only finite integrations of joint PDFs. Therefore, while a generic 
closed-form expression is not possible, the desired joint PDF can be easily numerically evaluated 
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with the help of integral tables [18], [19] or using standard mathematical packages such as 

Mathematica or Matlab etc. 

Theorem 5.3: (Joint PDF of j n]K and Yl In-.K) 

n=l n=m+l 

m K a 

Let X = In-.K and Y — In-.K, then we can simply obtain the joint PDF of Z = [X, Y] 

n=l n=m+l 

as 

Pz(x,y)=p m k 3 (x,y) 

n=l n = m + l 



Pm-i k b -i (x - z 2 ,z 2 ,y - z 4l z 4 ) dz 2 dz 4: , 

™ „ =1 n=m+l 

m 

for x > — y. (45) 

K s ~ m 

Proof: Omitted. ■ 
Note again that only the finite integrations of joint PDFs are involved. 

VI. Closed-form expressions for exponential RV case 

The above novel generic results are quite general and apply to any RVs. We now illustrate in 
this section some results for the i.i.d. exponential RV special case, where the PDF and the CDF 
of 7 are given by 

P (7) = 4 cx P [ -Z J , for 7 > 0, (46) 



7 V 7 
and 

P (7) = 1 - cxp , for 7 > 0, (47) 

respectively, where 7 is the common average. Therefore, ( flOl ), (fT2l) and ([141) specialize to 
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[e (z a , -S i )]' 

[c(z a , Si)Y 
[H (z a , Zb, -St)] 1 



exp 


;-( 








771 


ft)™ 




rn 



y> (-1) J (m 



E 

3=0 



exp 



exp 



[exp (-z a Si)] J 



Si 



(48) 



(49) 



exp (- ((m - j) z a + jzb) 5,) 



(™ - j) \ ( 3 
— = z a exp --z b 

7 / V 7 



(* + *) 



(50) 



After substituting (1481) . (l49l) and (l50l) into the derived expressions of the joint PDF of partial 
sums of ordered statistics presented in the previous sections, it is easy to derive the following 
closed-form expressions for the PDFs by applying the classical inverse Laplace transform pair 
given in (|60l ) and the property given in (loTI ). While some of these results have been derived using 
the successive conditioning approach previously, we list them here for the sake of convenience 
and completeness. 



A' 



1) PDF Of J2 ln:K. 

n=l 



K-l 



PZ (zi) 



_£l 



(51) 



A 



2) Joint PDF of "/ m:K and J2 In-.K- 



71=1 



(A 



-iy.{K-2)\^ K ex P ( Zl ^ 22 ) 



K-l 



x E (-!)'( i)[*2-3Zi] n ~'U(z 2 -j Zl ) form = l, 



Pz (zi,z 2 ) = < 



A! 



(A-m)!(m-l)!(A-2)!7 A ' 
A-rn 



exp 



(52) 



X 



E (-l) J rj m )[*2-(m + j-l)*x] 



A-2 



J=0 

xJ7(z 2 - (m + j - 1) zx) 



for m > 2. 
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m K 

3) Joint PDF of £ j n:K and J2 In-.K- 

n=l n=m+l 



A! 

(A-?ri)!(A-m-l)!(m-l)!(m-2)!7 



7T exp 



7 



x / dj m:K [z\ - m-f m:K ] m 2 U(z 1 ~ m-jm-.K) 
o 

A — m . T 

X E (-I)' ( K 7 m ) [«8 - hm:K] m U (z 2 - hm:K) 
3=0 



Pz (zi,z 2 ) = < 



F exp(-^) 



(A'-l)!(A'-2)! 7 -f 
J=0 



X 



A'., 



4) PDF Of E ln:K. 

n=l 



A! ( x \ 

(A-A 3 )!(A s -l)!(A s -2)!7»s CX P ^ — f J 



Pz> (x) = < 



xJo* 



l-exp(-^ 



A-A s 



TO > 



§ exp , 

7 ^ V 7 



where 



1 — exp — 



1 - exp ( -| 



m m 



A-l 

a: s = i 



= £(-ir 



i=0 









(7) 


exp ^ 


-;)] 









A'., 



5) Joint PDF of ^ m . K and E 7n:A-: 
a. For m = 1, 



71=1 



pz(zi,z 2 ,z 3 ) 



K\ 



(K-K s )\{K s -2)\{K s -m K > 

K~K„ 



exp 



Zl + Z 2 + ^3 

7 



Z3 

1 — exp — — 
' 7 



W(«i-z 3 ) 



A a -2 

><E 

J=0 



(-1) 

x U (z 2 - - 2 - j) z 3 - jzj) 
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b. For 1 < m < K s — 1, 



Pz (z ll z 2 ,z 3 ,z 4: ) 



K\ 



(K - K s )\ {K s -m- 1)! (K s - m - 2)! (m - 1)! (m - 2)!^ 



x exp 



2_1 + z 2 + z 3 + z 4 

7 



z 4 

1 — exp — — 

' 7 



[z\ - (m - 1) 2 2 f 



A' s — m— 1 



, /if s - m - 1\ .. . ,/c TO _2 

) [z 3 - (K s -m-l- j) Zi - jz 2 \ 



xU (z 2 - z^ U (zi - (m - 1) z 2 ) U (z 3 - (K s - m - 1 - j) z 4 - jz 2 ) 



c. For m — K s — 1, 



PZ (zi, 2 2 , 23 ) 



A"! 



(Jf - AT S )! (AT 5 - 2)! (AT S - 3)!t*- 



exp 



Zl + z 2 + z 3 



d. For m = K s , 

Pz{zi,z 2 ) = 



Z3 

1 — exp — — 

7 



xU( Zl -(K s -2)z 2 ). 



K\ 

(K - K s )\{K a -l)\(K s -2)\ri K ' 

K~K 3 



U (2 2 - z 3 ) [ Zl - (K s - 2) z 2 



K.-3 



exp 



zi 

1 — exp — — 
' 7 



zi + z 2 

7 



z 2 - (if, - 1) z^" 2 17 (z 2 - (K s - 1) Zi) 



(57) 



(58) 



(59) 



VII. Applications 

The above derived joint PDFs of partial sums of ordered statistics can be applied to the 
performance analysis of various wireless communication systems. In this section, we discuss 
two examples. 

A. Example 1) 

In conventional parallel multiuser scheduling schemes, a particular scheduled user's signal 
is detected by correlating signals of all scheduled users at the receiver. Therefore, under these 
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practical conditions, every scheduled user is interfering with every other scheduled user. These 
effect is called Multiple Access Interference (MAI) or Multiple User Interference (MUI). This 
MUI is a factor which limits the capacity and performance of multiuser systems and the level 
of interference becomes substantial as the number of scheduled users increases. To take into 
account the effect of MUI caused by the other scheduled users, the signal to interference plus 
noise ratio (SINR) which measures the ratio between the useful power and the amount of noise 
and interference generated by all the other scheduled users is used. The level of interference 
becomes substantial as the number of the scheduled users increases because the SINR can 
decrease considerably in these conditions. This decrease in SINR can lead to a certain reduction 
of the rate allocated to each scheduled user. With the above motivation in mind, the impact 
of interference on the performance (throughput) of the scheduled users assuming a selection 
based parallel multiuser scheduling scheme is needed to investigate the total average sum rate 
capacity and the average spectral efficiency (ASE) based on the SINR of the scheduled users. 
The major difficulty in investigating this total average sum rate capacity and ASE resides in the 
determination of the statistics of the SINR of the m-th scheduled user. Based on our approach, 
we can derive these results and then be applied to obtain the total average sum rate capacity 
and the ASE. Based on this system model, the statistics of the SINR of the m-th scheduled user 
can be derived using the joint PDF of SNR of m-th desired scheduled user and the sum of the 
SNRs of the interfering (K s — 1) scheduled users among total K (K > Kg) users. Application 
to our proposed approach and the correspondent performance results over was presented in [20]. 

B. Example 2) 

Minimum selection generalized selection combining (MS-GSC) is an adaptive diversity com- 
bining scheme [5], [21], [22]. The basic idea is to combine a minimum number of best diversity 
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paths among L ones such that the combiner output SNR is above a certain preselected threshold, 
denoted by jt- The performance analysis of MS-GSC scheme is challenging due to the ordering 
operation on diversity paths and adaptive combining operation. In [22], the average symbol error 
rate (SER) of MS-GSC was calculated as the weighted sum of the conditional average SER given 
that m paths are combined, with the weights being the probabilities of combining m paths and the 
average SER of conventional generalized selection combining (GSC), which always combines m 
best paths, used as the conditional average SER. However, because the receiver with MS-GSC 
may combine m best paths only under the condition that the combined SNR of first m — 1 best 
paths is below the output threshold, the statistics of combined SNR with MS-GSC given that m 
paths are combined is different from that with corresponding conventional GSC. As such, the 
average SER result in [22] should serve as an approximation. To obtain the exact statistics of 
the combiner output SNR with MS-GSC, we need the joint statistics of the mfh largest RV and 
the partial sum of the m — 1 largest RVs [5]. Specifically, the exact CDF of the combiner output 
SNR with MS-GSC involves the probability Pi[J2f=i 7 i:L < 7 T & 7 T < J2f=i ly-L < x], which 
can be expressed in terms of the joint PDF of -j m: L and X)j=i lj-L- With our proposed analytical 
framework, this joint PDF can be obtained in a systematic fashion, with the generic expression 
for generalized fading environments given in Eq. (36). This new statistical result allows us 
accurately evaluate the performance of MS-GSC scheme over general fading channels. 

Additionally, this proposed method can be applied to the performance analysis of various 
wireless communication systems with diversity techniques such that the performance analysis of 
GSC RAKE receiver with self-interference and so on. 



August 19, 2010 



DRAFT 



S.S. NAM et at: AN UNIFIED FRAMEWORK ... ON ORDERED STATISTICS 22 

Appendices 
Appendix I 

Useful inverse Laplace transform pair and property 

It is easy to see from (f56i (1571) . (1581) and <|59]) that the derivations of the PDF from the MGF 
involve the classical inverse Laplace transform pair [18] 

^{(-l-)"}^^"^- *>0,n = l,2,3,..., (60) 

and the Laplace transform property [18] 

C^ 1 {e- as F(s)} = f(t-a)U(t-a), a > 0. (61) 

Appendix II 
Derivation of I m 

In this appendix, we derive Eq. (flOl) . 

Let us first consider the case m = K. Noting that p (jk-.k) exp (X^k-.k) = d (ik-.k,X), we 
can rewrite Ik as 

1K-\:K 1K — 1:K 

J d'jK-.KP {IK-.k) exp (X^K-.k) = J djK-.Kc' (■JK:K,X) (62) 


= c^ K .. K ,\)\l K - 1 - K 

= C (jK-l:K, A) ■ 

For the case of m = K — 1, after applying integration by part and (|62l ), we have 

JK-2-.K JK-I-.K 


djK-l:Kc' (jK-UK, A)c (jK-l-.K, A) 

7K-2:K 
2 1K-2-.K r 

= [c^k-i-.k,^] - I dr/K-i-.Kc('yK-i:K,X)c (jk-i-.k,*)- (63) 

o 
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After moving the integration part of the right hand side (RHS) to the left hand side (LHS) and 
some manipulation, we can show 

7K-2:K 

J d-y K -i:KC (jk-i-.k, A) c (jk-i-.k , A) = ^ [c iciK-i-.K, A)] 2 . (64) 



Using (|64l) in ( |63l ), ( 1631 ) can be re- written as 

7^-2:^ 7K-1:K 

y d-f K ^ 1:K p (ik-i-.k) cxp (A7k_i:k) y d'jK-.Kp {jk-.k) exp (Atr-^) = ^ [c (tr--**-, A)] 2 . (65) 



Similarly for the case of m = K — 2, with the help of integration by part, we can obtain the 
following 

JK-3.K 7K-2:K 7K — 1:K 

y djK-2-.KP (ik-2-.k ) exp {\^ K -2:K ) y d'jK-i-.KP (7k-uk) exp (A7if_i : if) y d-y K ; K p(-f K:K )exp(\"{ K .. K ) 
ooo 



7K-3:K 



djK-2:Kc' (jK-2:K, A) - [c (7fi:-2 : K, A)]' 



1 



7if-3:if 



= ^ [c(7if-2:K, A)] 3 - / djK-2:K [c{^K-2:K,\)Y C (^ K -2:K,\)- (66) 



2 



7if — 3:if 







2 / , 







After some manipulation and substitution, we have 

lK~3:K 7K— 2:K JK—l-.K 

y dj K -2:KP(jK~2:K)cx-p(\j K -2:K) J d'jK-hKP (lfK-l:K) exp (XJK-I-.k) J dj K -.KP {jK-.K ) Cxp {Xjk-.k) 


= ^[c(7K-3:K,A)] 3 . (67) 
This process can be generalized to arbitrary m, which leads to the result in Eq. (flOl) . 

Appendix III 
Derivation of I" b 

In this appendix, we show the derivation of Eq.(fl4l). 
Let 

I'a,b= y d-y b - 1: K P{-fb-l:K)cxp(X-f b ^ 1:K ) y d76-2:KP(7&-2:x)cxp(A7b-2:K) 

7b:K 7b — 1:K 

7a:A' 7a:R' 

x y dj b ^ 3:K p(j b - 3:K )cxp(\-f b - 3 ; K ) ■ ■ ■ y d7o+i:jirjj(7o+i:Jf)exp(A7 a +i:ji-). (68) 

7b-2:A 7a + 2:ftT 
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Using similar manipulations to the ones used in the previous Appendices UH we can write 

7a :K 

d^a+l-.KP {la+l:I<) CXp (A7a+l:if) = (1 { r Ja+2:K- l la:K, A) , (69) 

7a + 2: K 

la:K 7a:K 

d-Ja+2-.KP (la+2:K) exp (Xj a+ 2:K) J ^7o+l:JfP (la+l-.x) CXp (Xja+l-.K:) = ^ ]p (la+3:K, 7a: if, A)] 2 , 

7a + 3:fiT 7a + 2:/f 

(70) 

7t»:K 7a:K 1a:K 

dja+3-.KP {ja+3-.K ) exp (A7 a+3: _ft:) J d'f a +2:KP {la+2-.K ) exp (A7 a+2: x) J d'fa+i-.xP {la+i-.x) exp (A7 a +l :K r) 

7a + 4:if 7o+3:K 7a + 2:/C 

[M (7a+4:if,7a:if, A)] J . (71) 



3x2 



Using these results, I" b can be found in closed-form as 



7o:K 7a:A' 

I'a,b= J dr/b-l-K p (76-l:if ) exp (A76_l:Jf ) y d7 6 _2:ifP(76-2:if)cxp(A7;,_2:if) 

7b:K 7i>— 1:K 

Ja:K Ja:K 

X y d7f,-3:ifP (7&-3:if ) exp (A76_3 : k) ■ • ' y d7a+l:ifP (7o+l:if ) exp (A7o+l:K-) 

lb-1.K la + 2:K 

1 



( b _ Q _ !)| If 1 (76:if,7a:if, A)] 



(6-a-l) 



(72) 



Appendix IV 
Derivation of (1221) 



Starting with (1221) . by simply applying (flOl) . we can obtain the following result easily 

lm:K IK-UK 

d7™+i:ifP (7™+i:if ) exp (A 2 7„ i+ i:if ) • • • / d^KiKP (jk-.k) exp (A 2 7if:if) =— -7 [c(7 m:J f, A 2 )] (if_m) . 

(A — mj! 





By inserting (1731 ) into (12TT) . the MGF has the following form: 

OO 71: K 

MGF Z (Ai,A 2 ) = F fdf 1:K p (7i:if ) exp (A 2 7i : if )Jd^2-KP (72:if ) exp (A 2 72:if ) 



7tti-2:K 

X • • • X / d7m-l:ifP (7m-l:if ) exp (A 2 7m-l:if ) 



(73) 





7m-l:K 

X / ^7m:ifP (7m:if ) exp (Ai7 TO:J f ) J— rr [c (7m:if, A 2 )] (K_m) . (74) 

7 (A - m)! 



August 19, 2010 DRAFT 



S.S. NAM et at: AN UNIFIED FRAMEWORK ... ON ORDERED STATISTICS 25 

By applying the integral solution presented in ©, we can re-write d74l) as the following: 

oo 

MGF Z (Ai,A 2 ) = F d-y m:K p(^ m:K )exp(X 1 -/ m:K ) — [c (j m -.K, A 2 )] (K_m) 

J (K-my. 
o 

oo oo 
X /d7m-l:ifP (7m-l:if) exp (A 2 7 m -l:i<r) J rf7m-2:XP(7m-2:Ar)cxp(A 2 7 m _ 2:A ') 

7ro:K 7m-l:if 

oo 

x ■■• x Jd'ji-.KP (ii-.k) cxp (X 2 ^i:k)- (75) 

72: K 

By simply applying (fT2j) . we can obtain the following result easily 



drfm-i-.KP (Tm-i.-Jf) exp (A 2 7 m _i : ^) / dj m - 2: KP (l m -2:i<) exp (A 2 7 TO - 2: ^ ) 

7m:Jf 7m-l:ff 
OO 

x ■ • ■ x Jdr/uKP (1i-.k) exp (A 2 7i:k) 

72 ;K 

= ^ L T y[e(j m :K,\2)] {m - 1) (76) 



By inserting (1761) into (1751) . we can obtain the second order MGF of Z 1 = 7 m: ^ and Z 2 

K 

in-.K easily as the following: 

n=l 



MGF Z (Ai, A 2 )=— r-; — / o?7m:ifP(7m:K) e xp(Ai7 m: i f ) [c(7m:K, A 2 )] (K m) [e (7 m: K, A 2 )] (m :) 

(A — mj! (m — lj! ' 



(77) 



Appendix V 
Derivation of (1281) 

Starting with (|28l ), by simply applying (flOl ), we can obtain the following result easily 

7m:K 7K-1:K 

y d^m+uKP (jm+i-.K ) exp (A 2 7 ot +i:.r:) • • ■ y dj K:K p (jk-.k) exp (A^^jr ) 



[c(7 m :A-,A 2 )] (K - m) . (78) 



(A' - m)\ 
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By inserting (|78|) into (|26l) . the MGF has the following form: 

oo 

MGF Z (Ai,A 2 ) = Ffd^uKP (ji-.K ) exp (XiJuk ) 
o 

7m-l:K 

x--x d-y m -.KP (im-.K) exp (\i"/ m:K ) — [c(7 m .jf, A 2 )] (K ~ m) . (79) 

7 (it -to)! 

o 

By applying the integral solution presented in ©, we can re- write ( 1791 ) as the following: 

oo 

MGF Z (X 1 ,X 2 ) = F d% n . K p{j m . K )exp(X 1 j m . K )— [c (j m:K , X 2 )] {K ~ m) 

J {K-m)\ 
o 

oo oo 

x J d'jm-i-.KP (7m-i:i<r) exp (Ai7 m _i :R r) • • J dju K p (7i:/c) exp (Ai7i : k ). (80) 
By simply applying (fT2l) . we can obtain the following result easily 



djm-i-.KP {im-i-.n) exp (Ai7 TO _i : jf) ■ • ■ J d^ 1:K p (7i : k) exp (Ai7i : k) 

7m:K 72:A' 

= r ^-^[e(7m:K,X 1 )}^- 1) . (81) 
(m — 1)! 

m 

Substituting (|8TT) in (1801) , we can obtain the second order MGF of Z-y = Y In-.K an d Z 2 = 

71=1 

is: 

E 7n:if as 
n=m+l 



MGF Z (Ai, A 2 ) 
if! 

(JiT-m)!(m- 1)! 



djm-.KP (7m:if) exp (Ai7 m: Rr) [c (7 m: K, A 2 )] (K m) [e (7„ 1: k, Ai)] (m . (82) 







Appendix VI 

Derivation of the joint PDF of ^ m ., K and Y In-.K among K ordered RVs 

n=l 

K s 

In this Appendix, we derive the joint PDF of 7 m: x and Y In-.K among K ordered RVs by 

n=l 
n=irn 

considering four cases i) m = 1, ii) 1 < m < K s — 1, iii) m = K s — 1 and iv) m = K s separately. 

m-l K s -1 

Consider first the case ii), 1 < m < K s — 1. Let Z\ — Y In-.K, Z 2 = J m -.K, Z 3 = Y In-.K 

n=l n=m+l 
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and Z 4 = jKs-.K- The 4-dimensional MGF of Z = [Z 1 , Z 2 , Z 3 , Z 4 ] is given by the expectation 

MGF Z (Ai, A 2 , As, M)=E {exp {\ 1 Z 1 + X 2 Z 2 + X 3 Z 3 + \ 4 Z 4 )} 
= F drfx-.KP {ii-.k) exp (Xi^i-.k) • ■ • / d^ m -x, K p {j m -i :K ) exp (Xi-fm-i-.x) 



o o 

7m-l:K 

X J d'Jm-.KP {im-.K) exp {Mlm-.x) 


7m:K 7K S -2:K 

X y rf7m+l:ArP(7m+l:if)cxp(A 3 7 m+ i : A:) • ■ • J d^K,-l:KP{lK s -l:K)e^£>{MlK a -X:K) 


lK s -\-.K 

x / d'y Ks:K p("f K3 ; K )eyL-p(\4 : 'y Ks:K ) [c{^ Ks :k)] {K ~ Ks) ■ (83) 



o 

With the help of ©, (fT0|). (fl2l) and (O, we can easily obtain the 4-dimensional MGF of 

m-l K s -1 

Z\= In-.K, Z 2 = r ) w ,:K, Z 3 = J2 j n:K and Z 4 = -y Ks ; K as 

n=l n=m+l 

MGFz (Ai, A2, A3, A4) 

OO 



rf7K s :KP(7^ s :K)exp(A47K s :K) [c(7K s :K)] (K ^ 



(K a - m - 1)! (m - 1)! 

o 

OO 

X / d7m:KP(7m:if)exp(A 2 7 m: x) [e(7 m: K,Al)] (, ™ _1) [^(7K s :AT,7m:K:^3)] (Ks_Tn ^ 1) - (84) 



Having a MGF expression given in (1841) . we are now in the position to derive the 4-dimensional 

771— 1 — 1 

joint PDF of Zx = J2 In-.K, Z 2 = J m -.K, Z 3 = ]T j n:K and Z 4 = j Ks -.k- Letting Ai = -Si, 

7i=l 7i=m+l 

771—1 

A2 = —S2, A3 = —S3, and A4 = — S4 we can derive the 4-dimensional PDF of Z\ = Jn-.K, 

71=1 
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K s -X 

Z 2 = Jm-.K, Z 3 = j n:K and Z 4 = Jk s -.k by applying an inverse Laplace transform yielding 

n=m+l 

Pz (zi, z 2 ,z 3 , z 4 ) = £ Si 1 s 2j 5 3j 5 4 {MGF Z {-Si, -S 2 , -S 3 , S4)} 

{K~K 3 ) 



(K s -m- 

00 

X J dj m :K 

1K S :K 



!)i( m _ !)! y d lK 3 -.KP (ik s -.k) C s l {exp (-5 4 7x 8 :i<r)}) [c (7kvk )] 


P (l m :K) Cgl i CX P (-^Tm:^)} C-sl {l e (7m:Jf. -5l)] (m_1) } 



rP (22) P (-24) [c(z 4 )] (Jf ■ ftrs) C/ (z 2 - Z4) 



-m - l)!(m - 1)! 

x£^ {[e (* a , -5 1 )] (m - 1) } £^ {[fx (z 4 , z 2 , -5 3 )] (/ ^-" i - 1) } . (85) 



With this 4-dimensional joint PDF, letting X = Z 2 and Y = Z\ + Z 3 + Z 4 we can obtain the 
2-dimensional joint PDF of Z' = [X, Y] by integrating over z\ and z 4 yielding 

P Z '(x,y)= / / pz (zi, x,y - z 4 ,z 4 ) dz x dz 4 , (86) 

JO J(m-l)x 

or equivalently we can obtain the 2-dimensional joint PDF of Z = [X, Y] by integrating over 
z 3 and z 4 giving 

px r{K s —m—X)x 

p z '(x,y)= / / p z (y - z 3 - z 4 ,x, z 3 ,z 4 ) dz 3 dz 4 . (87) 

JO J (K s -m-l)z4 

K s -X 

We now consider the case i) for which m — 1. Let Z\ = 7i : k, Z 2 = ^in-.K and Z 3 = Jk s -.k 

n=2 

for convenience. For this case, the 3-dimensional MGF of Z = [Z 1 ,Z 2 ,Z 3 ] is given by the 
expectation 

MGF Z (Ax, A 2 , A 3 ) = E {exp (\ X Z X + X 2 Z 2 + X 3 Z 3 )} 

OO 

= F d^x-KP (iuk) exp (Ai7i :K r) 





71 : K 



J d r y 2 :KP(72:K)cxp(X 2 ^2:K) ■ ■ ■ J d^K.-x-.KV (^k.-i-.k) exp (A27X,-X:X) 



7K S -1:K 

x / ^7A' s :A'P(7A' s: A')exp(A37K s: A') [c(7K s :A')] ( ' ft:_ ' ft:s) - (88) 
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With the help of and (fT4l) . we can easily obtain the 3-dimensional MGF of Z\ = Ji-.k, 

Ks-l 

Z 2 = J2 In-.K and Z 3 = j Ks . K as 

n=2 

oo 

MGF Z (Ai,A 2 ,A 3 ) = F I d^ Ks , K p (jk.-.k) cxp (Xs^k b -.k) [c{^1k 3 :k)] {K ~ Ks) 



xj drfuKP (li-.K ) exp (Xi^i-.k) ^ 2 (7if s :if , 7i:if, A 2 )] (gs 2) . (89) 



Having an expression from the MGF as given in (1891) . we are now in the position to derive 

K s -1 

the 3-dimensional joint PDF of Z\ = Ji-.k, %2 = Yl In-.K and Z 3 = ^Ks-.k- Letting Ai = —Si, 

n=2 

X 2 = —S 2 and A 3 = — S3, we can derive the 3-dimensional joint PDF of Z x = 71.^, Z 2 = 
In-.K and Z 3 = ^Ks-.k by applying an inverse Laplace transform yielding 



n=2 



p Z (z 1 ,Z 2 ,Z 3 ) = C s *s,,Sa i MGF Z (Si, S 2 , S3)} 

OO 

= F J dj Ks :KP{jK s :K) £33 i ex P(S3lK s :K)} [c(1K 3 :k)] {K ~ Kb) 



oo 

x J dyi,KP(~tv.K)Csl {° X P (-^171^)} J^~Jy C sl {\m (ik s -.k,1v.k, -5 2 )] (As_2) } 

djK a :KP (1K s -.k)5 (Z 3 -"/K a :K) [c (^K.-.k)] ^ 



1K S :K 

OO 



(if. -2)! 

u 

00 

x J dj 1:K p (7l;g) £ (zi - 7l:Jf ) ^ 2)! ^ {[M (TA-.r-g.TitX?-^)]^'"^} 

JK 3 :K 

{R /_ 2) , P (*i) P (23) [c (z 3 )] ( " s) (*i ~ 23) { [M (*a ,«!,-&)] (Ks ^ 2) } • (90) 



With these 3-dimensional joint PDF, letting X — Z\ and Y = Z 2 + ^3, we can obtain the 
2-dimensional joint PDF of Z' = [X, Y] by integrating over z 2 yielding 

AK a -2)x 

Pz'( x >y)= Pz (x,Z2,y - z 2 )dz 2 . (91) 

J (^)y 

K s -2 

We now consider the case iii) for which m — K s — 1. Let Z\ = Y j n: K, Z 2 = Jk^-i-.k an d 

n=l 
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Z 3 = JKs-.k- The 3-dimensional MGF of Z = [Zi, Z 2 , Z 3 ] is given by 

MGF Z (Ax, A 2 , A 3 ) = E {exp (\ 1 Z 1 + X 2 Z 2 + X 3 Z 3 )} 

OO 7K S -3:K 

= F J d-f 1 ., K p(j 1 ; K )cX-p(\ 1 J 1:K ) ■ ■ ■ J dj Ks ^2:KP(lK s ^2:K)cxp(X 1 J Ks ^2:K) 




K s -2 



With the help of © and (fT2l) . we can easily obtain the 3-dimensional MGF of Z\ = 'jn-.K, 



n=l 



Z 2 = lK s -v.K and Z 3 = ^ Ks , K as 

OO 

MGF Z (Ai, A 2 ,A 3 ) = F I d^ Ka , K p(^ Ks , K )cyLv(\ 3 ^ Ka , K )[c(^ Ks -.K)] {K ~ Ks) 



x J d~iK s -i:KV (ik.-i-.k) exp (A 2 7if s _i:if) J^~Zq)\ I e (lK s -i-.K, Ai)] (A ' s 2) . (93) 
Having a MGF expression as given in (1931) , we are now in the position to derive the 3- 

K s -2 

dimensional joint PDF of Z x = J2 ^ n ., K , Z 2 = 7k s -i-.k and Z 3 = Jk s :K- Letting Ai = -S^ 

n=l 

K s -2 

X 2 = —S 2 and A3 = —S3 we can derive the 3-dimensional joint PDF of Z\ = j n: K, 



n=l 



Z 2 = Jk s -1:K and Z 3 = ^k s -.k by applying an inverse Laplace transform giving 

Pz 0i, z 2 , z 3 ) = Csts2,s 3 i MGF z (Si, S 2 , S 3 )} 



F 



■\(K-K t ) 



, K _ 2 ^, J d-y Ks :Kp(lK s :K) {exp(-S 3 J Ks :K)}[c(-fK B :K)) ( 


OO 

x J dr/ K .-l:KP(lK.-l:K) {exp (S 2 7k s -1-.k)} £5* { [e (lKs-l:K, Sx)]^ K '~ 2 ^ 

1K S :K 

(K S F - 2 )\ P{Z2)P{Z3) [c^)] {K ~ K3) U(z2-z 3 )C- s l {[e(z 2 ,-Si)} iK °- 2) } . (94) 



With these 3-dimensional joint PDF, letting X = Z 2 and Y = Z\ + Z 3 we can obtain the 
2-dimensional joint PDF of Z = [X, Y] by integrating over z 3 yielding 

Pz' ( x >v) = Pz (y ~ z 3 ,x,z 3 )dz 3 , (95) 



11 
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or equivalently we can obtain the 2-dimensional joint PDF of Z = [X, Y] by integrating over 
Z\ yielding 

Pz 1 ( x >v) = [ pz(zi,x,y-zi)dzi. (96) 

J(K s -2)x 

K s -l 

Finally, we now consider the case iv) for which m = K s . Let Z\ = ^k^.k and Z 2 = In-.K- 

n=l 

For this case, the 2-dimensional MGF of Z = [Z 1 , Z 2 ] is given by the expectation 

MGF Z (Ai, A 2 ) = £{cxp(AiZi + A 2 Z 2 )} 

= F J dji: K p (1i-.k) ex P (^27i:if) • * ' / d^K.-i-.KP {ik.-i-.k) exp {Mik.-i-.k) 



7K S -1:K 

x / d~<fK s :KP (iKs-.k) exp (XiIKs-.k) {c{ik b :k)] {K ~ Kb) ■ (97) 







With the help of © and (fT2l) . we can easily obtain the 2-dimensional MGF of Z\ = 1k s -.k and 



Ks-l 



Z 2 = E ln:K aS 



n=l 

OO 



MGFz (Ai, A 2)=7^^T; y d-iK.JcPilK.-Jc) exp (Ai7jr.:jr) [c(7^:A')] (K_Ks) [e(7A' a :K,A 2 )] (Xs - 1) . (98) 



Having an MGF expression as given in (1981) , we are now in the position to derive the 3- 

K s -1 

dimensional joint PDF of Z\ = 1k 3 -.k and Z 2 = In-.K- Letting Ai = —Si and A 2 = —S 2 



n=l 

K s -1 



we can derive the 2-dimensional joint PDF of Z\ = 1k s -.k and Z 2 = ^ i n , K by applying an 



n=l 



inverse Laplace transform yielding 



Pz ( Zl ,z 2 ) = £ s ls 2 {MGF Z (Si, -S 2 )} 



F 



(i<s i; 



dlK 3 -.KP (ik s -.k)£ Si {exp (-SiJk s -.k)} [c (ik 3 -.k)] 



^{[e^-.K-S^-^] 



(K-K.) 



-s 2 
F 



(K s - 1)! 



p ( Zl ) [c ( Zl )] {K - K ° ] C£ {[e (z u -5 2 )] ( ^- 1} } . (99) 
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(a) Example of original M-dimensional groups 
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(b) Example of substituted (M + 1) -dimensional groups 
Fig. 1. Examples for 3-dimensional joint PDF with non-split groups. 
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(b) Example of substituted split groups 
Fig. 2. Examples for 3-dimensional joint PDF with split groups. 
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Derive Target 
M-dimensional 
MGF 



Applying 
inverse L.T. 



Consider 
subsituted 
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Consider 
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based on 
substituted 
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Applying 
inverse L.T. 




Applying 
inverse L.T. 
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Transfer the derived joint PDF to the 
M-dimensional joint PDF 



Target M-dimensional PDF 



Fig. 3. Flow chart to obtain the desired A/-dimensional joint PDF. 
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(a) For the case of m = 1 
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(b) For the case of 1 < m < K s — 1 
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(c) For the case of m = A" s — 1 
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(d) For the case of m = K a 



Fig. 4. Joint MGF of -y m:K and ^ 7«:if • 

n=l 



August 19, 2010 



DRAFT 



S.S. NAM et al: AN UNIFIED FRAMEWORK ... ON ORDERED STATISTICS 



37 



Sung Sik Nam (S'05, M'09) was born in Korea. He received the B.S. and M.S. degrees in Electronic Engineering from Hanyang 
University, Korea, in 1998 and 2000, respectively. Also He received the M.S. degree in Electrical Engineering from University of 
Southern California (USC), Los Angeles, CA, USA, in 2003, and the Ph.D. degree at Texas A&M University (TAMU), College 
Station, TX, USA, in 2009. From 1998 to 1999, he worked as a researcher at The Electronics & Telecommunication Research 
Institute (ETRI), Daejeon, Korea. From 2003 through 2004, he worked as a manager at the Korea Telecom Corporation (KT), 
New Business Planning Office, New Business Planning Team, Korea. Since June 2009, he has been with the Department of 
Electronic Engineering at Hanyang University, Korea. His research interests include the design and performance analysis of 
wireless communications, diversity techniques, power control, multiuser scheduling. 



Mohamed-Slim Alouini (S'94, M'98, SM'03, F'09) was born in Tunis, Tunisia. He received the Ph.D. degree in electrical 
engineering from the California Institute of Technology (Caltech), Pasadena, CA, USA, in 1998. He was with the department 
of Electrical and Computer Engineering of the University of Minnesota,Minneapolis, MN, USA, then with the Electrical and 
Computer Engineering Program at the Texas A&M University at Qatar, Education City, Doha, Qatar. Since June 2009, he has 
been a Professor of Electrical Engineering in the Division of Physical Sciences and Engineering at KAUST, Saudi Arabia., 
where his current research interests include the design and performance analysis of wireless communication systems. 

August 19, 2010 DRAFT 



S.S. NAM et al: AN UNIFIED FRAMEWORK ... ON ORDERED STATISTICS 



38 



Hong-Chuan Yang (S'OO, M'03, SM'07) was born in Liaoning, China. He received the B.E. degree from the Changchun 
Institute of Posts & Telecomm. (now part of Jilin University), Changchun, China, in 1995. During his graduate study at the 
University of Minnesota, USA, Dr. Yang received an M.Sc. degree in Applied and Computational Mathematics in 2000, an 
M.Sc. degree in Electrical Engineering in 2001, and a Ph.D. degree in Electrical Engineering in 2003. Dr. Yang has been with 
the Department of Electrical and Computer Engineering Department of the University of Victoria, Victoria, B.C., Canada since 
September 2003. He became an associate professor in July 2007. From 1995 to 1998, Dr. Yang was a Research Associate at 
the Science and Technology Information Center (STIC) of Ministry of Posts & Telecomm. (MPT), Beijing, China. His research 
work focuses on different aspects of wireless communications, with special emphasis on channel modeling, diversity techniques, 
system performance evaluation, cross-layer design and energy efficient communications. Dr. Yang is a recipient of the Doctoral 
Dissertation Fellowship (DDF) Award from the graduate school of the University of Minnesota for the 2002-2003 academic year. 
He is serving as an editor for IEEE Transactions on Wireless Communications and Journal on Communications and Networking. 
Dr. Yang is a registered professional engineer in B. C, Canada. 



August 19, 2010 



DRAFT 



